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Abstract 

There are many fields where the transition from diff'usive to ballistic motion is important. Here we deal with relaxation 
' processes in nmr in gases. Correlation functions for trajectory variables (position and velocity) valid across this transition are 

^—i known for several geometries in the case of specular wall reflections. In this work we show that the conditional probability 
, density p (1^, t|~j^o, io) for a random walk satisfies the telegrapher's equation and present an analytic solution for this function. 

^SJ ■ We will use it for calculating the relaxation due to an axion mediated force and a magnetic dipole impurity. 
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I. INTRODUCTION 

The influence of magnetic field inhomogeneities on nmr measurements has been discussed by many authors. It is 
usually sufficient to treat the problem using second order perturbation theory, (Redfield theory, [l]) according to 
which frequency shifts and relaxation rates can be shown to depend on the frequency spectrum of the magnetic field 
fiuctuations seen by the ensemble of particles moving through the inhomogeneous field. This spectrum is given by 
the Fourier transform of the correlation functions of the various field components. 

Usually the field variations are assumed to be described by constant gradients so the correlation functions of the 
field components are proportional to the correlation functions of the particle positions. McGregor 0] has shown how 
to calculate the correlation functions using diffusion theory. Recently in [3| it has been shown that these trajectory 
correlation functions can also be generalized for conditions (long mean free path for gas collisions) where diffusion 
theory is not valid and a general form has been given valid for all values of mean free path. 

The correlation functions are calculated by averaging the two successive positions, "af (t) , x' {t') over a joint proba- 
bility distribution which in turn is calculated in terms of the conditional probability distribution, p ^"a^, t\x' , t'^ which 

can be calculated by diffusion theory when that is valid 0] . By making use of the form of the conditional probability 
distribution calculated by diffusion theory it has been possible to apply these ideas to fields with arbitrary spatial 
variation [5] and the technique has been applied to placing a limit on a possible axion mediated unconventional force 
which falls off exponentially with distance from the walls of the container 

In this work we extend the work on the generalized correlation functions to calculate a generalized joint probability 
distribution valid for all values of mean free path. The generalized conditional probability is shown to satisfy a partial 
differential equation. This allows the study of the transition between diffusion and ballistic motion for the case of 
arbitrary field variation and should be applicable to a wide range of physical problems. 

In [3] the authors show that the velocity autocorrelation of confined gasses including gas collisions (wall collisions are 
assumed specular) can be expressed as a weighted sum of solutions to the stochastic damped harmonic oscillator and 
appropriate boundary conditions. They solved for the velocity correlation function of gasses confined in cylindrical 
vessels , and used its spectrum to predict a frequency shift important in searches for particle electric dipole moments 
, @ • In the authors applied the same idea to the rectangular cell combined with the observation that specular 
refiections enable a correlation in 2 dimensions to be decoupled into two correlations of 1 dimension. In 10] the 
authors derive the position autocorrelation function starting with the previously determined velocity autocorrelation. 
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II. FROM BALLISTIC TO DIFFUSIVE MOTION 



A. An equation for the conditional probability density 

A free (collision free) particle starting at xq with velocity v has a conditional probability density given by 

I 4,^0 = 0) = (5(5t + ^) (1) 
Equation ([IJ satisfies the wave equation as it is a function of {x — vt) . In three dimensions 

^ D V p[t,t I XQ.to) = (2) 

To include the effects of collisions we introduce a damping term as in [4] and write 

0^ + - Q-^ V V p{¥ ,t \ xo.tn) = Q (3) 

This is known as the Telegrapher's equation. It has applications to the propagation of sound in an absorbing 
medium, to a damped vibrating string, among others, and provides a method for eliminating the infinite propagation 
velocity from the theory of heat conduction or diffusion. 

For large ([3]) goes over to the wave equation, for small it becomes the diffusion equation. This allows us to 
study the transition region for intermediate values of Tc. Similarly for short times when the variation is rapid the 
second time derivative dominates showing that the short time motions is ballistic going over to diffusive motion for 
longer times when Tc is short enough. 

There is an extensive literature on the applications of this equation. We can only give some entry points to this. 
In reference the authors have taken a similar approach to ours in order to calculate the Green's function in 
a bounded region. The equation has been applied, among others, to the conduction of heat, since it overcomes the 
problem of the infinite propagation velocity of the Fourier theory [l^ and has been called the relativistic theory of 
heat conduction [l^. It is claimed to be applicable to heat conduction for distances on the order of or less than the 
collision mean free path [l^l, but this is contradicted by [l^. Chen's argument is specific to heat conduction and 
does not apply to the random walk problem considered here. A review of the applications of different forms of the 
equation to heat conduction is given in [l6|. Goldstein [l3], as early as 1951, in a little noticed paper, has discussed 
the Telegrapher's equation in relation to the random walk problem. 

It has also been applied to the transmission of light through opaque materials [l^ , [l^ . 

Here we apply it to the generic random walk problem and show how it can be used to study relaxation and frequency 
shifts in nmr. We show by comparing to simulations that our results apply for short distances and times 

B. A generalized conditional probability density 

Equation ([3]) can be solved by the usual separation of variables. Writing the space-time dependent function as 



F{t,t)/^-^ (4) 



we have 



The solution is a linear combination of terms with time dependence given by 

exp (-771,2*) (6) 

where 

m,2 = ^{^±s) (7) 

s = yjl - Auj'^t'^ (8) 

uj = kv (9) 
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We require a solution of ([5]) that satisfies the boundary conditions 

it ■^p{l^,t\^),to)sn^O, (10) 
yields the known diffusion theory result for p(l^,t \ xo,to) in the limit of small Tc, 

and yields a correlation function in agreement with the results of ref. 10], equation (26), where we now specialize 
to one dimension for simplicity. 

' (11) 

(12) 

(13) 



(r) - (a^ • 2^0)) , 

x{t) ■ xop{'^,t I xo,to)pixo,tQ) 



n=l,3,5. 



m 



4L2 



t ^ — ^ 1 



, , cosh - — t H snih - — t 

4 1 ^ 2t, / s,, \ 2r, 



(14) 



where 



s„ = Vl -4a;2r2 
1 

mrv 



m,2 



(l±s„) 



(15) 
(16) 

(17) 



and p{xo,tQ) is the position probability distribution at time t = tQ. 

Notice, for the ballistic limit, lim (s„) — ioo. and lim ( ) = icOn while for small Tc, Sr, 



1 — 2ujIt^ and 



2^ 



1 

2^ " 
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L J TD 



(18) 
(19) 



where tu is the time to diffuse a distance L. 

In addition p(x, t | a;o;0) rnust satisfy the initial conditions 

p{x, t = 0| xq, 0) = (5 (a; — xq) 
= 



dp{x,t\ 2:0,0) 



at 



(20) 
(21) 



t=o 



Strictly speaking the conditional probability distribution used here and in Q which satisfies the homogeneous 
equation with an initial condition equal to 5 il^) is not exactly the same as the Green's function which is the solution 
for a point impulse source. In the case of a linear equation containing only a first order time derivative it can be 
shown that the Green's function does have the required initial condition so that in the case of diffusion the functions 
are identical. This does not hold in the present case. 

The solution to ([3]) satisfying these criterion (in one dimension, — L/2 < a: < L/2) is 



1 I ■r-^ /"mTX\ /mTXct\ • fnnx\ . (n'KXo \ 1 

p{x,t\xQ,Q)^ -i 2^ cos j cos ^ j + 2^ sm j sm ^ j S 

yn=even n=odd ) 



X .. 



X < cosh 



Sjit 

2^ 



1 . , (Snt 

— smh I - — 



(22) 



We can see that this is a solution by considering the trigonometric functions in exponential form. Then every term 
will be of the form 



exp 



Tinx 
±1— ^- ?7i,2i 



(23) 
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so the equation is satisfied. The boundary and initial conditions are satisfied. The other requirements are shown to 
be satisfied by direct calculation. We take the sums over ± integers so as to include both directions of motion. 

The initial condition p{x,to\xo,to) = S{x — xq) is seen to hold since the inclusion of the terms multiplied by j- 
does not effect the initial condition, they are identically zero at t = 0. 

This solution applies to particles moving with a single velocity. It is necessary to average over the appropriate 
velocity distribution. The periodicity of the solutions means that reflections from the walls are properly taken into 
consideration (see appendix). 

This solution is seen to be correct in the limits of large and small Tc- The generalization to 2 and 3 dimensions is 
obvious. 

Correlation functions of arbitrary functions of position are seen to reduce immediately to combinations of the 
Fourier transforms of those arbitrary functions. 

In order to check its validity in the intermediate region we have performed Monte-Carlo simulations. The results 
are compared with equation (j22|) and the agreement is seen to be excellent. 



III. COMPARISON WITH MONTE CARLO SIMULATIONS 

We introduce dimensionless variables 



x'^^ (24) 
i' = - (25) 

Tc 



Then equation ^ can be written 



d'pi^^n Zt'o) ^ dp0,nZt'o) _ il^2^ij^^t I 5^,io) = (26) 

where Ic = vTc/ L — \c/ L is the normalized collision mean free path. In these units the velocity v' — Ic- We consider 
the solutions in the region —0.5 < x' < 0.5 and drop the primes from now on. 

The conditional density is simulated in one dimension. This is accomplished by tracking an ensemble of trajectories 
, A'^ = 500, 000, for a given initial position xq and velocity ±w. The trajectories scatter at random according to the 
cumulative distribution function of a collision in time, 

P — 1 — p"*/"^^ 

^ scatter ^ ^ 
^scatter '^c ln(l Pscatter^ 

where the probability of scattering, Pscatter, is a uniform distribution between 0, 1. At the start of each trajectory, or 
directly after a scattering event, a time of next scattering is selected by randomly selecting Pscatter- At the time of a 
scattering event there is a 50/50 chance of it being scattered in the reverse direction with speed v, versus maintaining 
its current trajectory. This is consistent with isotropic scattering in one dimension. Wall reflections are treated 
as specular. The whole treatment is easily modified to allow for the transport mean free path with non-isotropic 
scattering in higher dimensions. 

Below is a plot (fig.l) for = 0.05. This is in the region where diffusion theory is valid, however unlike the present 
theory the diffusion theory does not account for the peaks representing the unscattered particles moving with finite 
velocity. 

The Gibbs phenomena is clearly seen in this plot. The Gibbs phenomenon arises due to the impossible task of 
representing an abrupt change with a finite series of continuous functions. In the finite sum the abrupt change can 
only be approximated by the fastest oscillating component, and in the limit of an infinite series the fastest frequency 
is infinite and the exact form is recovered. However for purposes of plotting we cannot sum to infinite terms, and 
must accept the rapid oscillations, or smooth the rapid oscillations for a more pleasant and informative picture. To 
smooth the conditional density we take a running 4 point average. This is done for 4 x 500 points across the region. 
Typically 20,001 terms are summed prior to smoothing the functions but in this figure (fig. 2) we show the sum over 
only 301 terms for visibility. 

The smoothing correctly shows the average of the rapid oscillations and gives a clearer picture. The following 
(fig. 3) is a plot of the conditional density with = 0.3. This is in the transition region. The peaks representing the 
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Conditional Probability Showing Wall Reflection r=.05 



Monte Carlo t'=9 
Theory t'=9, 301 Terms 
Theory t'=1 2, 301 Terms 
Monte Carlo t'=12 




-0.5 -0.4 -0.3 -0.2 -0.1 



0.1 0.2 0.3 0.4 0.5 



FIG. 1: Conditional probability density with a;o = ,for the case Ic = 0.05, well into the diflusion limit, showing the agreement 
between the analytic result and the simulations for t' — 9 and 12. The peaks representing the unscattered particles are clearly 
visible at these times as is the Gibbs phenomenon due to cutting off the high frequency terms. 



Smoothed versus Unsmoolhed Conditional Density, i'^1 




FIG. 2: Conditional probability density with a;o = ,for the case Ic — l,t' = 1.1, approaching the ballistic limit. The particles 
have made one wall collision and have passed through each other {v' = Ic) . Shown are a truncated sum (A'^ = 301 terms) and 
the result of smoothing a sum with N=20,001. 



initial trajectory of the particles are seen to be large for longer than a collision time. It is seen that the theory agrees 
extremely well with the Monte- Carlo simulations. 

In the earlier time the trajectories are moving away from each other, in the later time they have reflected against 
the wall and are moving toward each other. The initial trajectories are scattering at a slow enough rate to leave a 
uniform wake of probability behind them, the uniform probability is deposited in layers behind the traversing peaks. 
That is to say the rate of change of probability going from the initial trajectory into the wake is slow compared to the 
rate of crossing the cell. The wake appears uniform and the curve reminiscent of diffusion is not seen. The following 
(fig.4)is a plot for /' = 1. More toward the ballistic region, however still somewhat in the transition region. 
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Smoothed 



Conditional Density r=0.3, 



N=20,000, 



4 - 
3 - 



- Monte Carlo t'=0.4 

- Theory t'=0.4 

- Sim t'=2.4 

- Theory t'=2.4 



FIG. 3: Conditional probability density with xq = ,for the case Ic = .3, t' = 0.4 and 2.4. A wall collision has taken place 
between the two times shown. The build up of the wake of scattered particles is clearly visible. 
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Smoothed Conditional Density versus IVIonte Carlo; 1=1 
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FIG. 4: Conditional probability density with a:o = ,for the case Ic — l,t' = 0.4, 3.2 and 100. By the later time the density is 
uniform and the unscattered peaks are gone, 



This behavior is similar to the I' — 0.3 case. The wake of probability deposited behind the peaks is less, however 
after a long time the peaks can no longer be seen and a uniform distribution is found. The proposed theory is shown 
to agree very well with the Monte- Carlo simulations. 

We see in general that as we are considering particles starting with velocities ±u, there are 2 peaks, representing 
unscattered particles moving away from xq in opposite directions. The peaks are decaying in time while a 'wake' of 
scattered particles builds up in the region behind them. The time development can be followed as the particles make 
multiple collisions with the walls. For short the wake morphs into the diffusion theory result, where the 'wake' is 
position dependent, after some time. For the ballistic case of longer Tc, the wake does not develop any significant 
spatial dependence while the peaks decay with time. In all cases the wake behind the unscattered particles are larger 
than the wake ahead of them. 
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IV. CALCULATION OF SPECTRUM OF CORRELATION FUNCTIONS. 



Now, starting with (1221) we can calculate the position correlation function according to (|T2)) and then take the 
Fourier transform to obtain the s pec trum. We can differentiate p2l) with respect to r to get the velocity position 
correlation function according to, [20| 

R.V (t) = (27) 

and a further differentiation will give the velocity autocorrelation function, [3| ■ We can calculate correlation functions 
involving quantities which vary arbitrarily with position. The results will be valid for short times and small distances 
allowing to consider fields varying over distances small compared to the collisional mean free path. 

As it is usually the spectrum of the correlation function that is of interest we give here the Fourier transform of 
that can then be integrated over arbitrary field distributions to give the desired spectrum. 



Gix,Xo,Lo)^ p{x,T \ XQ,0)e''^''dT (28) 
Jq 

4 { E (^) - m - E - (=F) - i"^) I ^ 



n — odd 



h (c.) ( 1 + + Qi {uj)(\--\) (30) 



Sn / \ S 



where 

Ql,2 W) = 2 , 2 31) 
^ + ^71,2 

are the Fourier transforms of e~^^-^'^ , with 771,2 given by ([TS]). We see the results for the correlation functions will be 
given in terms of Fourier transforms of the arbitrary field variations, similar to the case with diffusion theory [5], [6)]. 

This will now be applied to the problem of an axion mediated interaction with the walls, @ and to the problem of 
relaxation caused by a magnetic dipole impurity on the walls of the measuring cell, (2l| . 



V. DERIVATION OF THE FALSE ELECTRIC DIPOLE MOMENT SYSTEMATIC ERROR IN A LINEAR 
FIELD WITH THE GENERALIZED CONDITIONAL DENSITY. 

In the case of a constant gradient across the bound region, the proposed conditional density should give the same 
phase shift as |9'|. Using the equation (|27l) 

d 

{K {t - t) v4t)) = — {K{t)x{t + r)) 

Where the correlation of the field strength and position in terms of the conditional density for a linear gradient, 

— ¥ 

{hy{t)y(t + t)) = / dxoGxXo p{xot) / x p{x,t \ xo,t + t) dx (32) 



Where Gx is the strength of the linear gradient. 



n=odd \ c / n \ c 



{K (t-T) Vx{t)) 



SL^Gx 



E 



e sinh 



odd 



SnJ_ 

. 2t, . 



(34) 
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We now write it in a form which only sums n from 1 to oo.(we multiply by 2.) 



hxv{T) = (/la; - t) t):j(t)} = 



1 f s t\ 

^5 — sinh ( 7^ ) 



n—odd 



To find the phase shift we take the cosine transform according to the prescription in (Q) and 



Hxv (l^) — 







We notice that this field strength correlation spectrum gives the expected zero frequency value. 

TT* ^ — ' 

n—odd 

We arrive at the expected phase shift found in (3). 

duj{uJo) = — t;—Hxv{ujo) 
2c 
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VII. APPENDIX, CONSIDERATIONS OF WALL REFLECTIONS 



We set to = 0. This can be thought of as the time for the RF pulse, where the conditional density is a delta function 
and the initial probability density is flat with value . Now we consider a reflection from a single wall 

p{x, t\xo) =5 (^xO + vt-^ -{^ - x)j =S{xO + vt- L + x)) (35) 

The two terms should be summed for a complete description from t = 0. Now say we are confined from -L/2 .. 
L/2. So we make a second wall reflection giving us an impulse at 



and another wall 



and another. 



More Generally 
for V > and x > xo\to=o 



p{x,t\xo) = d{xO + vt-2L-x) (36) 



p{x,t \ xo) = 6{xo+vt -3L + x)) (37) 



p{x,t\xo) = d{xo+vt-4L-x) (38) 



p{x, t \ xo,to =0) = ^ S{xo +vt~ 2nL - x) + ^ 6{xo + vt- (2n + 1)L + x) (39) 



ra=0 n=0 

and for v > O\t=to and x < a;o|to=o 



p{x, t \ xo,0) = ^ 5{xo +vt- 2nL -x) + ^ 5{xo +vt- (2n + 1)L + x) (40) 



n=0 n=0 

and for v < O\t=to and x > Xo\to=o 



p{x, t\xo,Q) = 'Y^ 5{xq +vt + 2nL + a;) + ^ 5{xo +vt+ (2n + 1)L - x) (41) 
and for v < 0\t=to and x < xo\to=o 

oo oo 

p{x, t \ xo,0) = ^ S{xo +vt + 2nL + a;) + ^ 6{xo + vt+ (2n + 1)L - x). (42) 



n=0 n=0 



Now we make the assumption that at every xq there should be equal probability of a +v and —v. Furthermore 
we can combine the x < .To|to=o and x > .to|«o=o initial conditions. 

With that we have defined a conditional density for particles in one dimension with no gas scattering and purely 
specular wall reflections. We use the expansion of the delta function as: 

^(.-.0) = ^ E -K^)-^^) (43) 



+ 



n—even 
n=odd 



Where n goes from — oo to oo with even/odd respectively for each sum. 
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